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Abstract 

I discuss the effect of resummation of soft gluons in hadronic production of high 
mass systems, and in particular in heavy flavour production. I show that in widely 
used x-space resummation formulae, spurious terms that grow factorially in the 
order of the perturbative expansion are present. These terms spoil the convergence 
of the perturbative expansion. I also show that it is possible to perform the soft 
gluon resummation in such a way that these terms are not present. Implications for 
top and high mass dijet production at the Tevatron are discussed. 
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1 Introduction 



In this talk I will deal with the problem of the resummation of logarithmically enhanced 
effects in the vicinity of the threshold region in hard hadroproduction processes. Drell-Yan 
lepton pair production has been in the past the best studied example of this sort |1|, 0, ^ . 
The threshold region is reached when the invariant mass of the lepton pair approaches 
the total available energy. A large amount of theoretical and phenomenological work has 
been done on this subject. References @] and summarize all the theoretical progress 
performed in this field. Resummation formulae have also been used in estimating heavy 
flavour production [Q, [0. In this case only a leading logarithmic resummation formula 
is known. Calculations of the next-to-leading logarithms are in progress 0. 

Here I will not deal with sophisticated higher order effects calculations. I will in- 
stead describe some recent progress ^ in the understanding of how to implement the 
resummation. 



Coloured particles emit soft gluons with high probability. Normally the effect of soft 
gluon emission is small (at least in inclusive quantities), since they only slightly affect 
the kinematic of a process. However, in a production process of high mass objects, 
when we approach the threshold, soft gluon emission becomes important. Let us fix 
our attention on heavy flavour production in hadronic collisions near threshold. The 
process is schematically depicted in fig. ^ The incoming protons make a big effort in 



providing partons with a large fraction of the longitudinal momentum, thus going towards 
the very large x region of the structure functions. Under these circumstances, even the 



2 What are soft gluon effects 




Figure 1: Heavy flavour production near threshold 
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small amount of energy wasted by soft gluon radiation yields an important suppression 
of the cross section. In the usual application of the factorization theorem, part of the 
radiative corrections due to gluon radiation are included in the structure functions. Thus, 
depending upon the factorization scheme, the left-over suppression may be positive or 
negative. In the MS and DIS scheme the left-over is negative, so that the suppression 
effect appears instead as an enhancement of the cross section. 

Let us now focus upon the case of heavy flavour production. The perturbative expan- 
sion for the partonic cross section at order 0{a^) , neglecting obvious indices (incoming 
parton types, etc.), has the structure 



1 + Cas \og\l -p) + 0{as log(l - p), «f) 

where s is the square of the partonic centre-of-mass energy, p = Am? / s 
Born cross section 



(1) 

and cro(s) is the 
(2) 



The scale at which as is evaluated is as = as(m^), unless the argument is explicitly given. 
The reader can find in ref. explicit formulae for the functions f^^\p), as well as for 
the constant coefficient C. The precise value of C depends upon the type of incoming 
partons (i.e. quarks or gluons) and upon the factorization scheme. In both the MS and 
the DIS scheme it is positive, so that in the following we can think of it as being a positive 
constant. Resummation, according to ref. 0, gives 

„2 



^/W(p)exp 



Cas 



P) 



(3) 



where s' is a scheme-dependent function of s that goes to zero as s^Am?. We have 
s' = (1 — p)'^rn?, where r/ = 1 in the DIS scheme, and 3/2 in the MS scheme. Formula 
(^ is supposed to include all terms of order a^log"'(l — p) with n > m. Remember in 
fact that 

^s{s') = , / = — — ^2 ."^ , — -. 17 = tts (1 - asboT] log(l - p) + ...) . (4) 

Oology feolog^ + 6o^log(l -p) 

so that in the exponent there are terms with arbitrary powers of a^, and a power of the 
logarithm which is always larger than the power of as- The first subleading terms have 
the form a^ log''(l — p). 

In order to get a physical cross section the partonic cross section given above should 
be convoluted with parton luminosities 

^ J C{r)f'\p) exp [Cas ((1 - pfm^) \og\l - p)] dr 

(5) 
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where p = 4m^/(rS'), and (omitting obvious parton indices) 

C{t) = J F{xi)F{x2)S{xiX2 — t) dxidx2 ■ (6) 

Here we can spot a problem in the resummation formula. When performing our integral 
over r, as r — s> p the argument of in the exponential approaches zero. Before it actually 
hits the zero, it will hit a singularity in the running coupling as{s'), causing the integral 
to diverge. In order to avoid the divergence a cutoff fio was introduced in the literature 
|T^. Observe that this cutoff has nothing to do with the standard factorization and 
renormalization scale fi. It is essentially a cutoff on soft gluon radiation, imposed in order 
to avoid the blowing up of the running coupling associated with soft gluon emission. 

The use of a cutoff seems an ad hoc procedure in this case. It can however be justified 
to some extent. Suppose, for example, that we end up in a QCD calculation with a 
formula like 

G= / asik'^)Gik^)dk (7) 
Jo 

where G{k'^) is a smooth function as k'^—^0. Integrals of this kind are often found, for 
example, in the computation of shape variables in jet physics. This expression also is 
divergent as k'^ — >0, since at some point as approaches the Landau pole. The divergence 
can be handled by a cutoff /io, which has to be large enough for to be barely pertur- 
bative. For example, we may choose fiQ = 5A, a value around 2 GeV. We can then argue 
that ^ 

G= dk + G^. (8) 

In fact, the divergence of the 1-loop expression of as does not signal a real physical 
divergence. More likely, the point at which as becomes of order 1 signals the breakdown 
of perturbation theory. We therefore exclude this region, estimate its contribution by 
dimensional analysis, and obtain a power correction. A slightly more formal justification 
makes use of the concept of IR (infrared) renormalons. We expand eq. (|^) in powers of 
asiQ'^), using 
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where t = logQ^/fc^. The integral can be performed, and one gets 

oo 

G = 2asY.j\{2a,boy , (11) 

which is a divergent series, since a factorial grows faster than any power. This lack of 
convergence is in fact a general feature of the perturbative expansion in field theory. 
The perturbative expansion should be interpred as an asymptotic one. The terms of the 
expansion (|TT]) decrease for moderate values of j. As j grows the factorial takes over, the 
terms stop decreasing and begin to increase. This happens at the value of j at which the 
next term is equal to the current one 

(j + l)!(2aA)^+'=j!(2aA)^' (12) 

or roughly 

^ 2^0 ■ ^^^^ 
Asymptotic expansions are usually handled by summing their terms as long as they de- 
crease. Of course, in this resummation prescription there is an ambiguity, which is of the 
order of the size of the first neglected term. In our case 

j^J.{2asboy-- ^ e^i°s^-^l = ^ ^ , (14) 

f Q 

which gives a power-suppressed ambiguity with the same power law that we found using 
the cutoff procedure. 

From the discussion given above, we would expect that the resummation formulae 
should include a cutoff of the order of a typical hadronic scale, and varying the cutoff 
within a factor of order 1 should affect the cross section by terms of the order A/Q. In 
fact, this is not the case. The cutoff has a dramatic effect on the cross section, as can 
be seen from figs. 2 and 3 of ref. 0. For example, the uncertainty band obtained by 
varying the scale /io between 0.2m and 0.3m for top production in the gg channel brings 
about a change in the cross section of a factor of 2, for a top mass between 100 and 
200 GeV. These two values of fiQ correspond to cutting off the gluon radiation at energies 
of the order of 20 to 30 GeV, therefore much larger than a typical hadronic scale. Other 
proposals for the resummation procedure have appeared in the literature. In ref. a 



method was developed in the context of Drell-Yan pair production, and it was applied to 
the heavy flavour case in ref. [0. Also in ref. (as can be seen from formula (116) and 
the subsequent discussion) unphysically large cutoffs are present, much larger than the 
typical hadronic scale that one expects. 
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In the following section, I will show that the presence of large cutoffs and of large 
ambiguities in the resummation formula, is not at all related to the blowing up of the 
coupling constant. In other words, there are other sources of factorial growth of the 
perturbative expansion for the resummation of soft gluons, and they largely dominate 
the factorial growth due to the running coupling. I will also show that these terms are 
spurious, and that soft gluon resummation can be easily formulated in such a way that 
these terms are not present. 



3 Problems with the x space resummation formula 



For definiteness, let us focus upon the resummation formula (^. We pointed out earlier 
that this formula is divergent when the argument of as becomes too small, and the 
coupling constant blows up. In fact, formula (||) is divergent even for fixed coupling 
constant. At fixed coupling it can be written as 

„2 



^(-) = ^ /dr/:(r)/W(p)exp a,Clog2(l-p) 



a, 

m 



dpi ip/p) /W(p) exp kClog2(l - p) 



where p = Arn^/S, and the integral diverges as p— i>l, since the exponential 



exp(alog^(l -p)) = (1 -p) 



-alog(l-p) 



(15) 



(16) 



grows faster than any inverse power of 1 — p as p — > 1. This divergence can again be related 



to factorial growth in the perturbative expansion. Expanding formula (|T^) in powers of 
as we get 

„2 oo 



(res) 



a: 



oo 1 /-I/O 

E TiiCas)' / dp P-^C (p/p) /(°)(p) log^^l - p) 



(17) 



It is easy now to see that, due to its singularity for p— >1 the integral of log (1 — p) grows 
like (2/c)! . Let us make here the simplifying assumption that f^^\p) ~ 9{1 — p) ■ In a 
neighbour of the singularity, the integral has the form 



dp log 



2k I 



P) 



dt e~H 



t,2k 



'1 — e Jlogl/e 

where we have performed the substitution t = log 1/(1 — p) . The above integral equals 



/log 1/e 



dt e~H 



t,2k 



dt e~H 



t,2k 



log 1/e 



dte-H 



t,2k 



{2k)\ 



log 1/e 



dt e~H 



t,2k 



(19) 
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Since 



dt e-h^'' < (log l/ef^+^ (20) 



we see that the contribution to the integral near the singularity is dominated by the term 
{2k) \ . The power expansion for the cross section is then 

^(res) ^ ^ ^Ca,)' ^ ^ A:! {ACa,)' . (21) 

k=0 k=0 

The above formula is in fact appropriate for the case of heavy flavour cross section with 
a lower cut on the invariant mass of the pairQ. As in the previous example, the factorial 
growth of formula (21) will give rise to ambiguitites in the resummation of the perturbative 
expansion. These ambiguities are not, however, related to renormalons, since they occur 
also at fixed coupling constant. In fact, they are in general fractional powers of A/Q, where 
Q is the scale involved in the problem, and A is a typical hadronic scale. For example, in 
the case of heavy flavour production at fixed invariant mass of the heavy quark pair via 
the gluon-fusion mechanism, the ambiguity would have the form (A/Q)°'^^, and it would 
thus be extremely relevant even for very massive heavy flavour pairs. 

These large terms in the perturbative expansion are in fact spurious. They are an 
artefact of the x space resummation procedure. This can be easily understood with the 
following argument. Exponentiation of the gluon emission is possible because, roughly 
speaking, each soft gluon is emitted independently. This independence is however only 
approximate: the total momentum must be conserved. Momentum conservation, however, 
is a sub leading effect in the soft resummation formula. Yet, its violation leads to factorially 
growing terms. These terms are subleading from the point of view of the logarithmic 
behaviour, but very important from the point of view of the factorial growth of the 
perturbative expansion. The presence of large factorial terms due to momentum non- 
conservation can be understood also by simple arguments. The emission of k gluons, 
where each gluon has a limit on its energy Ei < i], leads to a phase space that is larger 
by a factor of k\ than the case when the total energy of emission is bounded J^Ei < t]. 
Thus phase space alone provides a k\ factor. We can see in more detail the origin of the 
{2k) ! term by considering the formula for the partonic cross section with two emitted soft 
gluons, implementing momentum conservation. We have 

(P) = o(2<^«s) ^ / — z — P '{p}6{p- ziZ2p)dzidz2dp . 

Z Tfl J i — Z]_ ^ i — Z2 _^ 

(22) 



■^A more accurate analysis, using the known behaviour of f^^^ near threshold would have yielded an 
extra factor of 4/9 in front of C in eq. (|2l|). 
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The leading logarithmic term of the above integral is given by 



a^'\p) = ^{2Ca,r^\og\l - p) f^\p) + . . . , (23) 

where terms with less than 4 powers of logarithms are neglected. We now see that the 
integral of the leading logarithmic term of a^'^\p) in p has a large factor ^ 4! due to the 
integral of log''(l — p), while the integral of the full expression, eq. (p2D , gives 



/^pa(^)(p) = l(2Ca.)^4 
JO 2 m"^ 




log(l - z) 
1 - z 








'dp'f^'\p) =0, (24) 



due to the vanishing of the integrals of the soft emission factors with the + prescription. 
In general, we see that if we take generic moments of a'^^\p) (i.e. / p"^dp a^''\p)), the 
leading log expression grows like (2fc)! , while the full expression grows only geometrically 
with k. 

The criticism described so far applies to the calculations of soft gluon effects in heavy 
flavour production given in refs. and 0. As one may expect, the large factorial terms 
give rise to large corrections to the cross section. Since the term of the perturbative 
expansion grow strongly with the order, they also give large uncertainties. In ref. 
the presence of large uncertainties is in fact recognized. In ref. it is claimed that the 
uncertainties are small. Even there, however, an unphysically large cutoff is needed in 
order to make sense out of the resummation formulae. 

A second, more subtle criticism of the resummation formulae has to do with the 



presence of 1/Q corrections that arise from infrared renormalons. It was shown in ref. |]l4 
that soft gluon resummation does not yield the correct renormalon structure of the Drell- 
Yan cross section. This proof was given in a simplified framework in which the renormalon 
structure can be computed exactly. This result suggests the absence of 1/Q corrections 
in Drell-Yan cross sections, an issue that is still much debated in the literature 
We would like to remark, however, that even iil/Q corrections were present in Drell-Yan 
and heavy flavour production, they would only be of the order of 1% for top production 
at the Tevatron. 

It is possible to formulate the resummation of soft gluons in such a way that the kine- 
matic constraints are explicitly satisfied, and no factorial growth arises in the perturbative 
expansion. It is enough to formulate the resummation problem in the Mellin transform 
space[|. 

^In fact, resummation formulae are usually derived in Mellin space. 
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In ref. we specify a resummation prescription that uses the Melhn space formula. 
We demonstrate that with such prescription there are no factorially growing terms in the 
resummed perturbative expansion. Yet, this prescription consistently includes soft effects. 
We call this prescription "Minimal Prescription" (MP from now on), because it does not 
introduce large terms that are not justified by the soft gluon approximation. 



4 Phenomenological applications 

We have computed various heavy flavour production cross sections using our MP formula. 
Details are given in refs. |10[- Here I report the most salient results. 

The importance of the resummation effects is illustrated in figs. |^, ^ and ^ where we 
plot the quantities 

'^gg "^qq (5gg + Sqq , . 

A99) ' (qq) ' (99) ,(m) ■ 

"^NLO C'NLO "^NLO "t C'NLO 

Here 5 is equal to our MP-resummed hadronic cross section in which the terms of order 
a1 and have been subtracted, and cr(NLO) is the full hadronic NLO cross section. The 
results for h at the Tevatron can be easily inferred from fig. ^, since the qq component is 
negligible at Tevatron energies. 

For top production, we see that in most configurations of practical interest, the con- 
tribution of resummation is very small, being of the order of 1% at the Tevatron. A 
complete review of top quark production at the Tevatron, based upon these findings, has 



already been given in ref. ||I0[- We also observe that, for top production at the LHC, soft 
gluon resummation effects are negligible. Of course, in this last case, there are other cor- 
rections, not included here, that may need to be considered. Typically, since the values of 
X involved are small in this configuration, one may have to worry about the resummation 



of small-x logarithmic effects |19 



We see from the figures that in most experimental configurations of interest these 
effects are fully negligible. One noticeable exception is h production at HERAb, at = 
39.2, where we find a 12% increase in the cross section. This correction is however well 
below the uncertainty due to higher order radiative effects. For example, from the NLO 
calculation with the MRSA' |^ parton densities and = 4.75 GeV, we get cr^^ = 



10.45 ^Jg^nb, a range obtained by varying the renormalization and factorization scales 
from m;,/2 to 2mb. Thus the upper band is 80% higher than the central value, to be 
compared with a 10% increase from the resummation effects. This result is much less 



dramatic than the results of ref. 
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Figure 2: Contribution of gluon resummation at order a'^ and higher, relative to 
the NLO result, for the individual channels and for the total, for bottom production 
as a function of the CM energy in pp collisions. 



5 Jet Cross Sections 

The interest in the effects of resummation on the behaviour of jet cross sections at large 
energy is prompted by the discrepancy between the single-inclusive jet-p^ distribution at 
large p^, as measured by CDF ||2^, and the result of the NLO QCD predictions [^. 



For simplicity we will study the effects of soft gluon resummation on the invariant mass 
distribution of the jet pair, which is, from a theoretical point of view, very close to the 
Drell-Yan pair production. Observe that other distributions, such as the Pt of the jet, 
have a rather different structure from the point of view of soft gluon resummation. In 
fact, while the jet pair mass is only affected by the energy degradation due to initial state 
radiation, the Pt of the jet may also be affected by the transverse momentum generated 
by initial state radiation, and by the broadening of the jet due to final state radiation. 

A study of the jet pair mass distribution is not of purely academic interest, since also 
for this variable an analogous discrepancy between data and theory has been observed 
[|^. Studies of resummation effects in the inclusive Pt spectrum of jets are in progress 
(M. Greco and P. Chiappetta, private communication). 
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Figure 3: Contribution of gluon resummation at order and higher, relative 
to the NL result, for the individual subprocesses and for the total, as a function 
of the top mass in pp collisions at 1.8 TeV. 
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Figure 4: Contribution of gluon resummaMon at order and higher, relative 
to the NLO result, for the individual channels and for the total, as a function of 
the CM energy in pp collisions. 
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In fig. ^ we show the following quantities: 

^(3) ^(3) .(3_) ^(3) , ^(3) , ^{3_) 

"gg "qg "qq "gg ^ "qg ^ "qq /o^^n 

^(2) ' ^(2) ' ^(2) ' ^(2) ^^"^ 

where 5^^-* is equal to our MP resummed hadronic cross section in which the terms of 
order have been subtracted, and a*^^^ = + a^^ + cXql^ is the full hadronic LO cross 
section (of order a^). We use as a reference renormalization and factorization scale for 
our results fi = Mjj/2. Notice that for large invariant masses the effects of higher orders 
are large. To understand how much is due to the first order corrections (which are exactly 
calculable [^) and how much is due to corrections of order and higher, we show the 



following quantities in fig. || 

<5(4) ^(4) 5(4) , ^(4) , ^(4) 

"gg "qg "qq "gg ^ "gg ^ "qq /r,«x 

a(3) ' cr(3) ' ' ' ^ ' 

where 6^^^ is now equal to the MP resummed hadronic cross section with terms of order 
subtracted, and a^^^ is an approximation to the full NL cross section, summed over all 
subprocesses, obtained by truncating the resummation formula at order a^. This figure 
shows that indeed most of the large K factor is due to the pure NLO corrections, with 
the resummation of higher order soft gluon effects contributing only an additional 10% 
at dijet masses of the order of 1 TeV. 

These results should only be taken as an indication of the order of magnitude of the 
correction, since we have not included here a study of the resummation effects on the 
determination of the parton densities. From this preliminary study it seems however 
unlikely that the full 30-50% excess reported by CDF for jet p^'s in the range 300- 
450 GeV could be explained by resummation effects in the hard process. It is possible 
that the remaining excess be due to the poor knowledge of the gluon parton densities at 
large x, sua idea pursued by the CTEQ group [EB]. 
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Figure 5: Contribution of gluon resummation at order and higher, relative 
to the LO result, for the invariant mass distribution of jet pairs at the Tevatron. 
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Figure 6: Contribution of gluon resummation at order a.g and higher, relative 
to the truncated 0{a^) result, for the invariant mass distribution of jet pairs at 
the Tevatron. 
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